Broken Ergodicity in a Stochastic Model with Condensation 
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We introduce a variant of the asymmetric random average process with continuous state variables 
where the maximal transport is restricted by a cutoff. For periodic boundary conditions, we show 
the existence of a phase transition between a pure high flow phase and a mixed phase, whereby the 
latter consists of a homogeneous high flow and a condensed low flow substate without translation 
invariance. The finite system alternates between these substates which both have diverging lifetimes 
in the thermodynamic limit, so ergodicity is broken in the infinite system. However, the scaling 
behaviour of the lifetimes in dependence of the system size is different due to different underlying 
flipping mechanisms. 
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The asymmetric random average process (ARAP) [1,2] 
represents a simple fundamental model for numerous ap- 
plications, especially in the field of interdisciplinary ques- 
tions. It has a close relationship to the q model describing 
force distributions in granular materials [3], aggregation 
and fragmentation processes (e.g. mass transport prob- 
lems) [4] and traffic flow theory, in particular a continu- 
ous Nagel-Schreckenberg model [5]. 

The truncated ARAP (TARAP) studied in this let- 
ter, although equipped with rather simple homogeneous 
short-range interaction and periodic boundary condi- 
tions, features the surprising behaviour of broken ergod- 
icity, i.e. a state space decomposition into dynamically 
unconnected subsets. In addition, these substates are 
totally different: the system may reside either in a ho- 
mogeneous high flow phase or in a condensed low flow 
phase where a finite fraction of the total mass is located 
on a randomly chosen site. While examples of infinite 
aggregation [4] and ergodicity breaking [6,7] are known 
for ID nonequilibrium systems, the simultaneous occur- 
rence of these phenomena seems to be undiscussed and 
relevant for a lot of applications. 

The ARAP is defined on a ID periodic lattice with L 
sites. Each site i carries a non-negative continuous mass 
variable m, GMq . In every discrete time step t — *■ t+ 1 for 
each site a random number 6 [0, 1] is generated from 
a time-independent probability density function (pdf) <fi, 
sometimes called fraction density, that may depend on 
the actual configuration m = {mi, . . . ,mi,)- So the uni- 
versal form can be written as (f> — 0(ri , . . . , mi, . . .). The 
fraction r% determines the amount of mass r^m, trans- 
ported from site i to site i + 1. The transport is com- 
pletely asymmetric, i.e. no mass moves in the opposite 
direction i + 1 — > i, and we get 



(1 - ri)mi + r i ^im i ^ 1 



(1) 



These update rules correspond to a parallel dynamics. 
Due to the conservation of the total mass M = J^i m i 
the density p — ^ is fixed. We would like to add that this 
so called stick representation is equivalent to a particle 
picture [1,2] not used in this letter. 



The simplest version of the ARAP is obtained if we use 
the state independent uniform distribution defined by the 
fraction density 0=1. We will refer to this system as 
the free ARAP. Some results have been derived for this 
model so far whereby the exactness of product measure 
for L — > oo and its form [1-3] 



P{m) = Y\ P{mi) with P(mj) 
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is the most relevant for our work. Eq. (2) has been ob- 
tained in the context of the q model [3] that shares many 
properties with the ARAP. Further results for the ARAP 
can be found in [1,2,8-10]. 

In the ARAP the mass rirrii transferred from site i to 
site i+1 is in principle unbounded. This is different in the 
truncated ARAP that we introduce here. In the TARAP 
all transfers of masses riirii larger than a cutoff A > are 
rejected. We like to focus on the simplest model only, the 
truncated free ARAP. The corresponding fraction density 
is then given by <f> = \\i <j>{ri, mi) with 

<t>{n,mi) = [1 - R{mi)\ S{n) + 9 (i?K) - n ) , (3) 

where is the Heaviside step-function and 



R{rrii 



min 1 



(4) 



represents the maximum possible fraction. Note that <fi 
has become locally state dependent, i.e. the pdf depends 
on the mass mi at the corresponding site explicitly. 

Without loss of generality we set p = 1 for the rest of 
the letter because every TARAP defined by (p, A) can 
be mapped onto a (1, ^-system (Fig. 1). Furthermore 
we introduce the rescaled cutoff 



2L"5A 



(5) 



which ensures L-independence of the critical point. 

We begin by investigating the relation between the 
steady state current J, defined by the average mass trans- 
fer per site, and the cutoff parameter A (Fig. 1). For 



1 



A — > the flow vanishes because the transferred mass 
per site is always smaller than A while for A — > oo the 
system behaves like a free ARAP and the value of J tends 
to its maximum J max = 1. 
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FIG. 1. Current-cutoff diagram for system size L = 100 
and different densities p (| = o, l = o, 2 = □). 

If we study the process in more detail we discover 
the system to exist in two phases whereby the transi- 
tion point A c = 1 has been determined numerically and 
by analytical approximations. 

For rescaled cutoffs A > A c the steady state mass 
distribution is nearly identical to the one of the free 
ARAP and approximately given by (2) . The flow is inde- 
pendent of A and corresponds to the maximum current 
"^high = J max - Therefore we refer to this parameter range 
as the high flow phase. 

For rescaled cutoffs smaller than A c the steady state 
of the finite system is a composition of two different sub- 
states: the system can either exist in a high flow state 
(with properties as described above) or a state given by 
a macroscopic condensate, i.e. a finite fraction of the to- 
tal mass M resides on one randomly chosen site (even 
in the thermodynamic limit). So the mass aggregation 
is proportional to L, but not extensive in space. The re- 
maining mass is distributed equally with an algebraically 
decaying mass distribution. The current of the conden- 
sate state depends on A but does not exceed J^ x = \- 
Therefore we call this state low flow state and denote the 
parameter region A < A c as mixed phase. 

In a finite system the transition probabilities between 
low and high flow states are small but nonzero. The 
system switches between these states while evolving in 
time and an alternating current-time relation is obtained 
(Fig. 2). Note that the switching time between the two 
states is much smaller than their lifetimes. 

This alternating behaviour is very similar to systems 
with spontaneously broken symmetry [6,7]. However, in 
general the flipping is between states of broken symmetry. 
For example in [6] an ASEP with two particles (+) and 
(— ) is introduced and in regimes of spontaneously broken 



symmetry the system switches between states dominated 
by (+) or (— ) particles. In case of the TARAP the flip- 
ping occurs between the symmetric high flow state which 
is translation invariant and a low flow state in which this 
symmetry is broken. 
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FIG. 2. Current-time diagram of the mixed flow phase. Pa- 
rameters used: L — 100 and A = 0.846. Each current value is 
averaged over 10 3 time steps. 

In the thermodynamic limit the average lifetimes tjj 
and T£ of the high and low flow states diverge (Figs. 3 
and 4). This implies that the steady state in the mixed 
phase is not unique. Ergodicity is broken and the steady 
state can either be in the high flow state or the low flow 
state depending on the initial condition. However, the 
lifetimes of the substates do not scale equally with size 
L (see next paragraph). 
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FIG. 3. Linear-logarithmic plot of the lifetime th of the 
high flow state in dependence of the system size L for several 
rescaled cutoffs A = 0.76(c), 0.80(D), 0.84(0) and 0.88(A). 

Although Monte Carlo simulations are difficult in the 
mixed phase since the lifetimes of the substates are very 
large it is safe to say that both the average lifetimes of 
low and high flow states diverge in the thermodynamic 
limit. Fig. 3 indicates that th increases exponentially for 
all A < A c . Assuming th ~ exp(a(A)i) we obtain that 
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FIG. 4. Log-log plot of the lifetime tl of the low flow state 
in dependence of the system size L for the rescaled cutoffs A 
= 0.92(d) and 0.88(D). 



a is an increasing function of A. Note that in the pure 
high flow phase A > A c the relation th = oo holds for 
any L. The corresponding data for t£ in the low flow 
state are better fitted algebraically than exponentially 
(Fig. 4). We see again that the exponents obtained by 
the assumption tl ~ increase while moving into 

the mixed phase, i.e. for A tending to zero. 

Although both transitions are driven by fluctuations, 
the above scaling behaviours reflect different switching 
mechanisms. While the high — ► low transition is based 
on a collective effect which involves all sites of the lattice, 
the opposite transition involves only the lattice site where 
the macroscopic condensate is located. Under these as- 
sumptions we have analytically derived approximations 
that reflect qualitatively the lifetime behaviour [11]. 

Note that in spontaneous symmetry breaking both 
transitions obey the same purely fluctuation driven mech- 
anism yielding the same exponential behaviour of the 
lifetimes. A flipping mechanism between two alternating 
symmetric substates that is not purely driven by fluctua- 
tions is studied in [12], resulting in a logarithmic lifetime 
dependence. 

The above results are mainly based on Monte Carlo 
simulations. We continue now by presenting a few ana- 
lytical approximations which support the numerical find- 
ings. 

Starting point is a closer look at the average mass 
transfer per site (local current) 



J(m) = (rra) 



for sC m < A 
for A < m < oo 



(6) 



This shows that the average mass shift J{m) tends to zero 
for m — > oo and m — » 0. So high (low) columns shrink 
(grow) very slowly and accordingly low flow states, resp. 
one-stick configurations, are very stable. On the other 



hand, homogeneous configurations maximize the current. 
We will exemplify this in the following paragraphs. 

First we study the low flow state using the approxima- 
tion 



{J{m)) P =J{{m) P% ) 



(7) 



Here Pi denotes the steady state single-site distribution 
of site i. Equation (7) holds for stable distributions 
Pi{m) ~ S(m — nii) which we assume here. Introduc- 
ing the nomenclature mi = (m)p i we obtain from the 
continuity equation in the steady state 



^iow = J {mi) = J{m i+1 ) 



for all i. 



(8) 



Restricting now to the case where only one high col- 
umn exists at site j, i.e. m + = mj > A, the re- 
maining columns are all of the same mass (see (8)), i.e. 
to_ = rrii < A with i ^ j. So, using (6) and (8), the 
quantities m± are related by m + m_ = A 2 and the mass 
conservation law pL = {L — l)m_ + m + . 

Equipped with these relations we are able to compute 
m± and finally the current in the low flow state Ji ow - 
Assuming L » 1 our calculations lead to the formula 



rmax 

— J low 



{l - v 7 ! - A 2 } . 



(9) 



So the current of the low flow state increases from zero to 
J™^ x for 0<A^ A c . For rescaled cutoffs larger than A c 
relation (9) is not defined and the one-stick model is not 
valid anymore. Thus we have derived an upper bound for 
the occurrence of the low flow state in accordance with 
the simulations. As Fig. 5 shows the analytical result 
agrees very well with Monte Carlo data. 
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FIG. 5. Low flow state current in dependence of the 
rescaled cutoff. Comparison of analytical (-) and numerical 
(L = 100(c), 200(D), 500(O)) data. 

By extending the calculations above to N mass aggre- 
gates [11] the validity range in terms of A shrinks and 
we obtain smaller critical cutoffs A c . Furthermore the 
flux increases. This indicates that the mixed phase is 
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described best by the one-condensate picture, especially 
in the vicinity of the critical point. This is confirmed by 
numerical investigations that have shown that configura- 
tions with two or more aggregates are not stable. The 
pre-condensates struggle for masses until only one stick is 
left. This mechanism of competition is different to coars- 
ening where aggregates emerge at several sites and then 
coalesce [4]. 

Focussing on the high flow state, resp. phase, every 
site carries the same mass p on average. For large L 
we approximate, motivated by the numerical results, the 
mass distribution by the product measure solution of the 
free ARAP (2) and calculate by the help of (6) the high 
flow state current 

•/high = {l - (1 + i^A)e- iiA j ^ J££ (10) 

for all rescaled cutoffs. So the high flow state exists both 
in the mixed phase and in the high flow phase. Equa- 
tion (10) also predicts that the current in the high flow 
state differs from only for small absolute cutoffs A 
in finite systems. We have verified this by Monte Carlo 
simulations although it is no longer possible to distin- 
guish the states by their flows which are nearly identical 
in this parameter range. So we used the appearance of a 
macroscopic condensate as a criterion. 

A model that is strongly related to the TARAP is stud- 
ied in [5,14] where a similar kind of phase separation 
as in the mixed domain has been observed. The un- 
derlying process is the continuum limit (Krauss model) 
of the Nagel-Schreckenberg cellular automaton model 
[13] that is defined by continuous velocities and spa- 
tial coordinates. Although this traffic model is given 
by more complex dynamics than the TARAP, both pro- 
cesses have a common feature: moves may be rejected if 
a uniformly distributed random variable exceeds a given 
threshold. So the TARAP can be viewed as a toy model 
that catches some of the fundamental physics behind the 
Krauss model. We have also rewritten a simplified ver- 
sion of the Krauss model in terms of fraction densities 
to point out the similarities with the TARAP explicitly 
[11]. 

In the Krauss model an additional congested phase has 
been observed [5]. An analogous phase, corresponding 
to a pure condensed phase, could be expected for the 
TARAP in the limit of small cutoffs. However, in that 
regime the lifetime of the low flow state is much larger 
than the corresponding lifetime of the high flow state, in 
particular for small systems. Therefore, in Monte Carlo 
simulations it is difficult to distinguish between a pure 
condensed phase and the mixed regime. Furthermore, 
our analytical calculations have not shown any evidence 
for a transition point separating a (new) pure low flow 
phase and the mixed phase. 

What are the essential ingredients of the TARAP lead- 
ing to the observed phenomena? First truncation and an 



unbounded local state space allow for an unlimited con- 
densation. The formation of infinite aggregates can also 
be found in a model of aggregation and fragmentation 
[4,15] or a special zero-range process [16], both equipped 
with unbounded but discrete state variables. However, in 
these examples the condensed phase is unique and not as- 
sociated with a high flow counterpart, i.e. high flow phase 
and congested phase are separated by a phase transition. 
The same is true for a discretized version of the TARAP 
that we have investigated [11]. 

Therefore we believe that in case of the TARAP the 
continuous state variables allow for the coexistence of 
high and low flow state in the mixed phase, yielding 
nonsymmetric ergodicity breaking in the thermodynamic 
limit. The TARAP is free of an intrinsic mass scale 
in contrast to models defined on integer state space, 
equipped with a smallest mass unit 1. Thus, in the 
TARAP we obtain a trivial dependence on the density 
p which is also reflected in the scaling law (p, A) <-> 
(l,A/p). However, in discrete systems physics may 
change in the low density regime, e.g. an ARAP with to- 
tal mass M = 1 is nothing else than a one-particle ASEP. 
This difference may indicate why in the continuous vari- 
ant the homogeneous high flow state can persist over the 
whole parameter range, resulting in the coexistence phe- 
nomenon. 

We would like to thank Joachim Krug, Kavita Jain and 
Martin Evans for interesting und helpful discussions. 
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